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IV. INTRODUCTION 


An aerodynamic phenomena which sometimes occurs is a variation of 
an aerodynamic force with control surface deflection which depends on control 
direction of motion of the surface. This effect, which can appear as 
illustrated in the following sketch, is referred to as "control 
hysteresis." 

Arrows indicate direction 
of control deflection. 



Using wind-tunnel data from early shuttle configurations and simplified wing- 
body shapes, J. Peter Reding and Lars E. Ericson of the Lockheed Missiles 
and Space Company, under a contract for the NASA Johnson Space Center, 
identified "control hysteresis" as a possible problem area for the space 
shuttle arbiter's entry (reference 1). They showed that regions of shock- 
induced separated flow might be expected to exist on the leeward wing 
surface, and that upward deflected control surfaces affect the extent of 
the separation. Bow shock-control surface shock interactions can take 
place ahead of tr.e elevon deflected downward into the wind. Also, 
movement of a control surface into a wing leading edge, wing fillet, or 
nose vortex can cause the vortex to burst. Lags in re-establishing flow 
field equilibrium may result from these phenomena and thus "control 
hysteresis" or a dependency upon the direction of control surface 
deflection may exist. Before discussing the possible ramifications of 
"control hysteresis," an examination of the orbiter's control 
philosophy is needed. 



In the normal operational mode, the space shuttle orbiter entry is 
directed by on-board computers from deorbit through landing. Attitude 
is control led\y both aerodynamic control surfaces and a reaction control 
system (RCS) consisting of twenty (20) 4000 N (900 lb) thrust, hypergolic-fueled 
rockets. The aerodynamic surfaces include elevens (acting as both 
elevators and ailerons), rudder, speed brakes, and body flap. The speed 
brake and body flap deflections are determined from preset tables. The 
elevens, as elevators, are used for longitudinal control; as ailerons, 
they are used either for turn coordination (angle of attack (a) > 18* or 
Mach number (M) >5), or for roll attitude control (a < 18“ and M < 5) . 

When the ailerons arc being used for roll attitude control, the rudder is 
used for turn coordination, otherwise it is inoperative. The RCS provides 
Tolling, pitching, and yawing moments. The roll RCS, operative until a 
dynamic pressure of 479 Pa (10 psf) is sensed, aids the ailerons in turn 
coordination. The pitch RCS, operative until a dynamic pressure of 
9S8 Pa (20 psf) is sensed, aids the elevators in longitudinal control. 

The yaw RCS is used in two modes. For a > 18* or M > 5, the yaw RCS 
is used for roll attitude control. This is done by using the yaw thrusters 
to produce sideslip and then allowing the orbiter' s positive dihedral 
effect to roll the vehicle. When a < 18* and M ^ S, the yaw RCS augments 
the rudder. If "control hysteresis" is present in the ailerons, the most 
active aerodyn^unic control during entry, two possible situations exist, 
depending on the flight conditions. If a > 18* or M > 5, the ailerons 
may be unable to properly coordinate the turn, because as the ailerons 
deflect, moving first one direction and then the other, "hysteresis" 
will result in different aerodynamic moments for the same deflection angle. 
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In addition, at these flight conditions, a lateral 
trim network is fed to the commanded aileron circuit to minimize RCS 
propellant consumption due to lateral center-of-gravity offsets and 
cross-wind effects. This network causes slight uncoordination of a roll 
maneuver resulting in the ailerons overshooting the proper position and 
therefore requiring a reversal in the direction of the surface movement. 
Thus, the error signal in the commanded aileron circuit will oscillate 
about zero. If tliis oscillation is large** than any built-in deadband, 
the ailerons will also oscillate. Since the turn coordination signal 
is also fed to the roll thrusters, while they are operational, and yaw 
rate and roll angle error signals are fed to the yaw thrusters, 
"hysteresis" in the ailerons aerodynamic moments could result in 
alternate (positive then negative) thruster firings, and sharply increase 
fuel consumption. iVhen the ailerons are used for roll control (a < 18^ 
and MSS), hysteresis could result in the orbiter continuously hunting 
to reach the commanded roll angle, i.e. orbiter oscillating about the 
commanded roll angle, and since turn coordination is the responsibility 
of the ludder augmented by the yaw thrusters, this oscillation cou d 
produce alternate yaw thruster firings. Thus, the major effect of 
"control hysteresis" was expected to be an increase in the amount of RCS 
fuel required for entry. Since the amount of fuel that is carried is 
limited, this increase could result in poor control of the orbiter, and 
limit its maneuverability. 

Wind-tunnel tests using a remotely controlled eleven model were 
conducted in the NASA- Langley Research Center's Continuous Flow Hypersonic 
Tunnel at a Mach number of 10.3. ;.,*ta, taken at 40 frames/sec as the 
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elevens were driven through their operational range, revealed some 
evidence of hysteresis. The low Strouhal number (the model elevon rate 
of seven (7) deg/sec corresponds to a full-scale rate of only 0.175 
deg/sec), and low Reynolds number (1.0 x 10^ based on length for the model 
as compared to approximately 5.0 x 10^ full scale), made the determination 
of the amount of full-scale "control hysteresis" impossible. 

Consequently, a six (6)-degree-of-freedom analysis was conducted to 
determine the system tolerance to hysteresis. 

Hysteresis was modeled by offsetting and/or modifying the slope of 
the roll due to aileron such that the values of the rolling moment changed 
with the direction of control surface travel. 

In mid 1973, development of a simulator to examine automated space- 
craft entries was begun by the author, et al, specifically to examine 
this and other problems associated with the space shuttle orbiter's 
entry. It is a six (6) -degree-of- freedom, intercctive simulator known 
as the Automatic Flight Dynamics Simulator, ARFDS (reference 2). \RFDS 
was modified to include "control hysteresis" and an analysis of the 
entry was made to examine the impact of this potential problem. 

This analysis is the subject of this thesis. 



V. SPACE SHUTTLE SYSTEM DESCRIPTION 


The physical characteristics of the space shuttle that were used 
are summarized in Table 1. The guidance scheme utilized is described 
in Appendix A, and the control system is described in Appendix B. The 
control and guidance scheme are applicable from deorbit to the Terminal 
Area Energy Management (TAEM) interface which occurs at 457.2 m/sec 
(1500 fps) and an altitude of 21.3 km (70 000 ft). 

The entry in the automatic mode is directed entirely by on-board 
computers. The guidance system software produces a series of angle-of- 
attack and roll attitude commands that the control system software utilizes 
to direct the RCS and surface deflections. 



VI. MIND TUNNEL RESULTS 


VindotonzMl tests ttere ccnducted in the 'IAS& Laogley Research Cent«r* 
Gantianous Flow Rjpersoaic Tvemel of a space shuttle oirt>iter aDdel with 
reatotel3r coatxolled ele'vans (both left and ele'vcns nove together). 

Rata were taken at UO Acsaes/sec. as the eleecn was dri'ven between its 
linits (-1^'* to 1 ^* to -1^**) at seieen (7) degs/sec. Figure 1 shows 
pitching Boaent (C^) data fhirings fron these tests tdiich indicate that 
"control hjsteresis" nej be present. The "control hjsteresis" observed 
in the data is very close to the balance accuracy; however, these trends 
repeated for all cases that were run, indicating a strong possibility 
that "control hysteresis" did exist. 

Te&t paraneters including a low Stroiihal nunber (sodel eleven rate 
scaled to a Ihll scale rate of only 0.175 deg/sec.) and a low Reynolds 
madber (1.0 x 10 based on length as coopared to $,0 x lO” fall scale) 
Bade it iapossible to predict the aaount of fhll -scale “control hyste- 
hysteresis" that might be ejqiected. Therefore, sone analysis of the 
system tolerance to control hysteresis was needed. 
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VII. AUTOMA.TIC REENTRY FLIQIT DYNAMICS 
SIMULATOR (MIFDS) DESCRIPTION 

In nid 1973, developnent was begun on a new interactive sisulator 
capable of handling studies of automated spacecraft entries, namely, 
space shuttle. Up to that time, studies had been made at the NASA 
Johnson Space Center (JSC) and by the author using a batch mode ccmputer 
program known as G.E. Mass.* The new simulator was developed by 
personnel of the Analysis and Computation Division of the NASA Langley 
Research Center after being provided suitable equations and overall 
program logic by the author. The simulator is knoim as the Automatic 
Reentry Flight Dynamics Simulator, ARFDS, (see reference 2), and it 
has the following advantages over G.E. Mass: 

1. Interactive capability. 

2. Paster integration and table look-up capability. 

3. Separate time channels for guidance, control, and the 
equations of >K>tion. 

4. Entry states can be observed on strip charts as well as 
printed output. 

The dieckout of ARFDS was the responsibility of the author. To do 
this, the G.E. Mass Program was modified to simulate a typical shuttle entry 
problem. The author independently programed the control scheme, adapted a 


*Nritten under contract to the Manned Spacecraft Center <^nd was 
available in 1963. 
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guidance system routine* obtained from the shuttle program office at JSC, added 
an appropriate viscous aerodynamic model, added a reaction control system model 
complete with aerodynamic interference effects, and modified the integration 
routines to allow for a one (1) pass method for integration of the filters in 
the control system. 

The analysis described herein was performed on the checked out 
version of ARFDS to take advantage of the speed, ease of modification, 
and strip charts. For this analysis, models to simulate aerodynamic 
control hysteresis were added to ARFDS. 


*The guidance routine was rewritten for ARFDS to increase efficiency - 
see Appendix A. 
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VIII. SIMULATION APPROACH 

Since nominal entry runs had indicated that the elevens are more 
active as ailerons than as elevens, an investigation was made with ARFDS where 
"control hysteresis" was added to the rolling moment (C^) due to aileron 
deflection (5^) • This moment was chosen because: 

(1) Possible causes of "control hysteresis" such as separation and 
vortex burst would predominantly effect wing lift or normal force vdiioh 
determines the amount of rolling moment produced. 

(2) The rolling moment is the primary control output of the aile3K>n 
for which the control system is designed. 

"Control hysteresis" was added to Cj^ due to in two (2) different 
manners as shown in figure 2. Figure 2a shows the first method where 
if the aileron deflection is increasing (6^ positive), the n ominal valtie 
of is used; if the deflection is decreasing (d^ negative), the . linal 

value of is multiplied by Cl. This changes the slope of the Cg vs. 

^ values. Cl was varied from 0.5 to 1.5* providing a 50 percent 
change in the slope. For small values of d^. Cl would have little effect 
on CffOB this method allows for no displacement of the origin. To correct 
this deficiency, an increment, C2 (figure 2b), was added to the nominal value 
if ^ was negative. C2 was e^^ressed as an equivalent aileron deflection, 
e.g. C2 was equal to x" of aileron and therefore C2 was a function of M and a. 

C2 was varied between + 5" of aileron (the maximum d experienced in any 
of these cases was less than 10"). Figure 2o shows the range of hysteresis 
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covered by this study for Mach four (4) and Mach numbers greater than 
ten (j J). This range was felt to be greater than would be experienced 
by rl’<. orbiter. 

The aileron deflection remains zero (0) until a dynamic pressure 
(q) of 96 Pa (2 psf) is reached (nominally 320 seconds after deoibit). 
The simulation on ARFDS was initiated at this point and continued until 
the " XEM interface conditions are encountered. 



IX. DISCUSSION OF RESULTS 


The results of simulations for the nominal conditions (Cl s 1» 

C2 = O) with a guidance system sampling time of O.32 seconds and a 
oontrol system sampling time of O.OU seconds is shown on time histoi^ 
strip charts shown in figure 3. As indicated previously and in Ajipendix 
the ailerons axe used for turn coordination until a is reduced to 
18*** ^17^ seconds from deorhit. After that time* the ailerons are used 
for 4 control* and the rudder assumes the coordination role. Times of 
increased aileron activity, namely between I4OO and $00 seconds w3«ere 4 
increases from -1$** to >75** and after 1500 seconds where roll reversals 
occur, are times of increased RCS firings. Thus, it was expected that 
adding hysteresis would increase the BCS firings at these points and re- 
sult in more fuel consunqption. Figures U and 5 show the effect of re- 
ducing Cl to 0.7 and 0.5* re^eotively. Neither of these cases are 
significantly different from the nominal, either in the 6^ or the RCS 
thruster history. Returning Cl to its nominal value of 1.0 and applying 
C2 values of -2" (figure S) and -5^ (figure 7) produced similar results 
as did the combination of Cl s 0.5 and C2 « -5 (figure 8). On a batoh 
version of ARFDS (no strip charts)* cases for Cl = 1.5 with C2 = 0* and 
Cl s 1 with C2 s 5* >>sre also run. Table II shows the RCS lUel consunqption 
for all these oases and they fall within six (6) percent of the nominal. 

During the time of the present study, it became evident to the NASA 
Johnson Space Center and the contractor that the on-board-computers could 
become over-burdened by the many assigned tasks (flight control, guidance 
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and navigation functions, redundancy management, etc.)* It appeared 
desirable therefore to require guidance calculations Infrequently as 
possible. Guidance system design studies at NASA Johnson Space Center 
had shown that guidance system sampling every 2.0 seconds was sufficient 
for proper targeting. When the guidance system sampling time was 
increased from 0.32 seconds to 2.0 seconds, the author, et al, in 
reference 3, found that the RCS yaw thrusters would limit cycle causing 
the RCS fuel consumption to double. Three (3) system changes were 
identified that eliminate this limit cycling and actually result in a 
37 percent reduction in fuel consunqption from 320 seconds in the entry 
to the TAB1 interface over the nominal case identified earlier. The 
three system changes are: 

1. Replacement of the step changes in commanded angle-of-attack 
and roll attitude with a linear variation (ramp- like). 

2. Modification of two (2) gains in the control circuit. 

Figures 9-14 show the results of applying "control hysteresis" to the 
simulation with the revised system, and again "control hysteresis" makes 
little difference in control surface and vehicle motion time histories. 
Table III shows the fuel consumption due to hysteresis with revised 
guidance. 

Examination of the time history strip charts for both the nominal 
(0.32 second san^ling time) and revised guidance (2.00 seconds sampling 
time with system change) cases reveals that for most of the entry, 6^ 
is zero (0). Examination of the block diagram, figure 15, reveals 

that this must be due to the error signal remaining within the deadband 
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for most of the entry, and this is probably preventing the "control 
hysteresis" from having any significant effect. Consequently, this 
filter was removed from the circuit, and both a nominal (Cl » 1, C2 <> 0) 
case, and a "control hysteresis" case (Cl ° 0.5, C2 « -S) were run 
(figures 16 and 17) with the nominal guidance (0.32 second sampling time). 
The aileron without the deadband becomes a continuously moving control, 
and hysteresis results in significant increases in RCS fuel consunqption 
(see Table IV, where a Cl of 0.5 and C2 of -5 results in a 37 percent 
increase in fuel consunption over the Cl of 1.0 and C2 of 0 case). 

Thus, the deadband filter in the aileron circuit is acting to suppress 
the effects of any "control hysteresis" over the range of Cl and C2 
tested. There is no deadband filter in the elevator circuit, but since 
the elevator moves so slowly and pitch rates are quite small during the 
nominal entry, there is probably no problem. If, however, due to mission 
change, or in the presence of wind gusts, etc., the elevator becomes a 
more active control, a deadband filter will probably be required to 
suppress the effects of control hysteresis. This possibility should 
be the subject of a future study. 
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X. CONCLUDING REMARKS 

A study by J. Peter Reding and Lars E. Ericson of the Lockheed 
Missiles and Space Company identified "control hysteresis" as a possible 
problem area for the space shuttle orbiter's entry. This "control 
hysteresis" phenomena was detected in the NASA Langley Research Center's 
Continuous Flow Hypersonic Tunnel using a remotely controlled eleven model. 
Consequently, a six (6) -degree-of- freedom- simulation examination was 
conducted on the space shuttle orbiter's entry to determine if "control 
hysteresis" in the rolling moment due to aileron deflection has any 
major effects on the controllability of the orbiter or required reaction 
control system fuel consumption. "Control hysteresis" was modeled by 
offsetting and/or modifying the slope of the roll due to aileron such that 
the value of the rolling moment changed with the direction of control 
surface travel. 

The simulations indicated that the orbiter system can tolerate 
"control hysteresis" that produces a 50 percent change in the nominal 
characteristics, or an offset in the nominal characteristics equivalent 
to a five (5) degree aileron deflection with little increase in the 
required reaction control system's fuel consumption. This tolerance was 
traced to a deadband filter in the commanded aileron circuit. Removal 
of this filter results in significant increases in fuel consumption. 

Thus, "control hysteresis" in roll due to aileron is not a problem with 
the present control system for the range of hysteresis magnitude tested 
If conditions exist that make the elevator a more active control, a dead- 
band filter should be added to suppress the effects of "control 
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hysteresis." A study should be made to determine if a deadband filter 
should be added to the elevator signal. 
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TABLE I.- PHYSICAL CHARACTERISTICS OF SPACE SHUTTLE ORBITER 


Mass Properties 


MASS 

83 001 kg (182 986 lb) 

^XX 

1 029 066 kg-m^ (759 000 slug-ft^) 

lyy 

7 816 290 kg-m^ (5 765 000 slug-ft^) 

izz 

8 015 596 kg-m^ (5 912 000 slug-ft^) 

^XZ 

177 612 kg-m^ (131 000 slug-ft^) 

IXY « Iy2 = 0 



WING 

Reference Area 

249.91 (2690.0 ft^) 

Chord 

12.06 m (39.57 ft) 

Span 

23.79 m (78.06 ft) 


ELEVON 

Reference Area 

19.51 (210.0 ft^) 

Chord 

2.30 m (7.56 ft) 


RUDDER 

Reference Area 

9.30 (100.15 ft^) 

Chord 

1.86 m (6.1 ft) 

• 

BODY FLAP 

Reference Area 

12.54 (135.0 ft^) 

Chord 

2.06 m (6.75 ft) 
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TABLE II 

RCS FUEL CONSUMPTION FOR THE NOMINAL GUIDANCE 
WITH A SAMPLING TIME OF 0.32 SECONDS 


Cl 

C2 

FUEL CONSUMPTION, kg (lb) 

1 

0 

176 

(388) 

0.7 

0 

171 

(378) 

O.S 

0 

185 

(408) 

1.5 

0 

174 

(385)* 

1 

-2 

173 

(382) 

1 

-5 

181 

(399) 

1 

♦5 

176 

(389)* 

O.S 

-5 

178 

(393) 


•Cases were run on batch version of ARFDS so strip charts 
are not available. 
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TABLE III 

RCS FUEL CONSUMPTION FOR THE REVISED GUIDANCE 
WITH A SAMPLING TIME OF 2.00 SECONDS 


Cl C2 FUEL CONSUMPTION, kg (lb) 


1 0 

0.7 0 

0.5 0 


1 -2 


1 -5 


111 (245) 

112 (246) 

110 (243) 

111 (245) 
127 (280) 
124 (274) 


0.5 


-5 
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TABLE IV 

BCS FUEL CONSUMPTION FOR THE NOMINAL GUIDANCE 
WITH A SAMPLING TIIC OF 0.32 SECONDS WITH THE 
AILERON DEADBAND FILTER REMOVED 


Cl C2 FUEL CONSUMPTION, kg (lb) 

10 1S6 (343) 

O.S -S 213 (470) 








Figure 2. Hysteresis Models. 


AREA BETWEEN DASHED LINES INDICATES RANGE DP 
**CONTROL HYSTERESIS" STUDIED. 



(c) Range of "control hysteresis" studied 
Figure 2.- Concluded. 
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XU. APPENDIX A 

Analytic Drag Control Entry Guidance System 

The baseline guidance scheme controls the entry by roll modulation 
while flying a preselected angle-of-attack profile. Downrange is 
controlled by the magnitude of the roll angle and crossrange is controlled 
by multiple bank reversals. The guidance system outputs to the control 
system are commanded roll angle and commanded angle of attack. 

The Analytic Drag Control Entry Guidance System (ADC, ref. 4) was 
developed by the NASA Johnson Space Center to approximate an optimum 
entry profile determined previously. This profile is achieved by 
dividing the entry into five major phases as illustrated in figure A-1: 

1. Constant attitude phase 

2. Constant heat rate phase 

3. Equilibrium glide phase 

4. Constant drag phase 

5. Transition phase 

The space shuttle orbiter is commanded to fly a constant attitude 
trajectory until a specified total acceleration is attained. At 
this point, a constant stagnation heat rate trajectory is flown 
through pullout to a relative velocity of 6248.4 ra/sec (20 500 fps) 
or until the reference drag level becomes larger than that required to 
reach the target. If the latter condition is reached, the guidance 
scheme jumps to the constant drag phase. If this condition is not met. 
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an equilibrium glide profile is flown until either it intersects the 
constant drag profile required to reach the target and jump to the constant 
drag phase, or the velocity drops off to 2743.2 o/sec (9000 fps) . Whenever 
the velocity drops to 2743.2 m/sec (9000 fps), the transition phase is 
entered during which the commanded angle of attack is decreased to the 
value required at the Terminal Area Energy Management (TAEM) point, which 
occurs at a velocity of 457.2 m/sec (1500 fps) and an altitude of approxi- 
mately 21 km (70 000 ft) . 

Table A-I shows the input constants that were used, and figure A-2 
shows the block diagram of the guidance laws of reference 2 as 
modified for the Automatic Reentry Flight Dynamics Simulator. 

SYMBOLS 


PARAMETER 

UNIT 

DEFINITION 

AK 

sec"^ 

dD/dV for constant heat rate phase, 
used to define C 3 

ALDREF 

n.d. 

(L/D)j,g£, used in controller 

ALFM 

m/sec^ 

reference equilibrium glide drag 

ALMNl 

rad 

minimum roll command outside of lateral 
deadband (YB) 

ALMN2 

rad 

minimum roll command inside of lateral 
deadband (YB) 

ALPCMD 

deg 

o^, angle of attack command 

ARC 

m 

distance from intersection with alinement 
circle to target 

ARC 

rad 

(l/D)y/(L/Dj , used in roll command 
equations 
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ATK 

m 

radius of earth 

EA 

rad 

equilibrium glide roll angle used in 
iteration loop 

BAD 

deg 

final equilibrium glide roll angle 

BAl 

deg 

first iteration equilibrium glide roll 
angle 

BA2 

deg 

second iteration equilibrium glide roll 
angle 

CAGI 

2 2 
sec /m 

temporary calculation used in transition 
phase to calculate ALDREF and RDTREF 

CIGAR 

n.d. 

transformation matrix from Earth Centered 
Inertial (ECI) axes to geocentric axes 

C0SBADD 

n.d. 

temporary calculation in equilibrium glide 
ranging phase used to calculate DREFP 

CTH 

rad 

great circle range from orbiter to target 

C4 

m/sec 

parameter used to calculate RDTREF 

C5 

n.d. 

parameter used to calculate RDTREF 

Cll 

m~^ 

parameter used to calculate RERl and RDTREF 

C16 

2 

sec /m 

parameter used to calculate L^Dl 

C17 

sec/m 

parameter used to calculate ElDDl 

C21 

m/sec^ 

parameter used to calculate DREFP, RDTREF, 
SQ and TTll 

C22 

sec"^ 

parameter used to calculate DREFP, El, E2, 
RDTREF, SQ, TTll. and TT22 

C23 


parameter used to calculate C22, DREFP, 
El. E2, SQ, TTll, and TT22 

D 

N 

total drag force 

DBAR 

m 

distance from runway to alineraent circle 

DBB 

deg 

increment in roll angle in equilibrium 
glide phase 
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DELAZ 

rad 

azimuth error 

DF 

2 

m/sec 

final drag level in transition phase 

OLIM 

2 

m/sec 

control system limit drag level in 
transition phase 

DRAG 

m/sec^ 

current drag acceleration level 

DREFP 

m/sec2 

drag reference used in controller 

DT 

m 

planar range to target 

DTO 

rad 

angle between alinement circle center and 
tangency point 

DTR 

rad/deg 

It/ 180 

DVHEAD 

rad 

azimuth between runway and heading to 
tangency point of alinement circle 

D23 

m/sec^ 

parameter used to calculate AK 

EEF 

2 2 
m /sec 

current energy level 

EEF4 

m^/sec^ 

reference energy level used in transition 
phase 

El 

n.d. 

parameter used to calculate TT22 

E2 

n.d. 

parameter used to calculate TT22 

GAMMA 

rad 

flight path angle 

GCLAT 

rad 

orbiter geocentric latitude 

GCLATT 

rad 

target geocentric latitude 

GP 

m/sec^ 

acceleration of gravity at sea level 

GSTART 

n.d. 

acceleration in "g's" required to 
initiate constant heat rate phase 

HA 

m 

current altitude 

HAogrr 

m/sec 

dHA/dt 

IffiSER 

m^/ sec 

parameter in oblate earth correction term 
to RDTREF 

HS 

m 

altitude scale height 
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IDFG2 

n.d. 

switching flag in constant drag phase 

IDFG3 

n.d. 

switching flag in transition phase 

IFT 

n.d. 

initialization flag in equilibrium glide 
phase 

ISLECT 

n.d. 

phase selector 

ISTART 

n.d. 

initialization flag 

ISTP 

n.d. 

iteration flag in equilibrium glide phase 

ISTRT 

n.d. 

flag indicating acceleration level equal 
to GSTART has been reached 

ITR 

n.d. 

iteration flag in transition phase 

L/D 

n.d. 

lift to drag ratio 

(L/D)v 

n.d. 

lift to drag ratio in vertical plane 


n.d. 

minimum value of L0D1 

L0O1 

n.d. 

desired (L/Dy) 

PS IE 

rad 

current heading of orbiter 

PSIET 

rad 

current heading to target 

RAZ 

rad 

runway azimuth 

RCG 

m 

predicated range in constant drag phase 

ROC 

n.d. 

parameter used in RDTREF calculation 

ROTOLD 

m/sec 

final RDTREF in equilibrium glide phase 

RDT0L2 

m/sec 

final RDTREF in constant drag phase 

RDTREF 

m/sec 

altitude rate reference 

REC 

n.d. 

vector defining rvjiway coordinate system 

RECl 

n.d. 

[REC]-1 

REH 

m 

distance from center of earth to vehicle 

REQ 

m 

predicted equilibrium glide phase range 

RERl 

m 

parameter in range prediction for transition 
phase 
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RFF 

m 

predicted range in constant heat rate phase 

RG 

m 

vector from orbiter to runway center 

RGP 

m 

vector from orbiter to alincment circle 
center 

RK2R9L 

n.d. 

roll direction (♦ right, - left) 

RL0N 

rad 

orbiter 's longitude 

RLONT 

rad 

target longitude 

ROLLC 

rad 

roll angle command 

RPT 

m 

desired range in transition phase 

RPTl 

m 

range bias below velocity of 456.2 m/sec 

RTE 

m 

radius of earth at runway 

RTURN 

m 

radius of alinement circle 

Rll 

m 

first iteration of range prediction in 
equilibrium glide and transition phases 

R12 

m 

second iteration of range prediction in 
equilibrium glide and transition phases 

SQ 

sec*^ 

parameter used in constant heat rate 
range prediction 

SQQ 

sec*^ 

parameter used in constant heat rate 
range prediction 

TA 

m 

vector from alinement circle tangency point 
to vehicle 

TAP 

m 

vector TA in geocentric coordinates 

TARE 

m 

target vector from alinement circle 
center to runway 

TDREF 

m/sec^ 

parameter used in DREFP calculation in 
equilibrium glide phase 

TEMP 

m 

temporary calculation in equilibrium glide 
phase 

TRANCE 

m 

great circle range from orbiter to target 



59 


TTll 

m 

parameter used in range prediction in 
constant heat rate phase 

TT22 

m 

parameter used in range prediction in 
constant heat rate phase 

li 

2 

m/sec 

parameter used in calculation of ALDREF 

T2 

2 

m/sec 

constant drag level required to reach 
the target 

U 

rad 

1 DVHEAD 1 

OTARE 

n.d. 

TARE unit vector 

UXYZE 

’^.d. 

RG unit vector 

V 

m/sec 

earth relative velocity 

VBB 

m/sec 

intersection velocity between constant 
heat rate phase and equilibrium glide phase 

VCG 

m/sec 

predicted intersection velocity between 
constant drag phase and equilibrium glide 



phase 

VINERT 

m/sec 

inertial velocity 

V0LD 

m/sec 

final velocity in equilibrium glide phase 

V0LD2 

m/sec 

final velocity in constant drag phase 

VQ 

m 

predicted final velocity for constant drag 
phase 

VSAT 

m/sec 

reference circular orbit velocity 

VSATS 

m^/sec^ 

(VSAT)^ 

V10LD 

m/sec 

value of V/LD-152.4 

V20LD 

m/sec 

value of V0LD2-152.4 

XLFAC 

m/sec^ 

total acceleration 

XL0D 

n.d. 

L/D 

XYZE 

m 

geocentric position vector 

YB 

rad 

lateral deadband (amount of overshoot that 
guidance system will allow before commanding 
roll reversal) 



TABLE A- I.- ADC GUIDANCE INPUT CONSTANTS 


PARAMETER 

VALUE 

UNIT 

ALFM 

7.62 

m/sec^ 

ALMNI 

0.7986355 

n.d. 

ALMN2 

0.9659258262 

n.d. 

ATK 

6.36670702 x 10^ 

m 

DEAR 

14360.4 

m 

DF 

5.819 

2 

m/sec 

EEF4 

1.8580608 X 10^ 

2 

m /sec 

GCLA'fT 

34.55577617 

deg 

GS 

9.815 

m/sec^ 

GSTART 

0.05 

n.d. 

RAZ 

-0.7679448709 

rad 

RL0NT 

-120.5338 

deg 

RPT 

4.218856 X 10^ 

ID 

RPTl 

23150 

m 

RTE 

6.373298953 x 10^ 

m 

RTURN 

4632.96 

m 

VSAT 

7853.54 

m/sec 

VQ 

2133.6 

m/sec 
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Figure A-1.- Analytic Drag Control Entry Cjuidance System flow diagram. 
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REC(l.l) 

REC(1,2) 

REC(1,3) 

REC(2,1) 

REC(2.2) 

REC(2.3) 

REC(3»1) 

REC(3,2) 

REC{3,3) 


} INITIALIZE 


DATA GS,ATK,GSTART,RPT, 
VSAT.VQ.ALMNl.ALFM, 
RPTl.RAZ.RTE.DBAR, 
RTURN,ALMN2.EEF4,DF. 
RL0NT.GCLATT 


-CeS(RAZ)*SIN(GCLATT)*C0S(RL0NT)-SlN(RAZ)*SIN(RL0NT) 

-C0S(RAZ)*SIN(GCLATT)*SIN(RL0NT)+SIN(RAZ)*C0S(RL0NT) 

C0S(RAZ)*C0S(GCLAn) 

S I N ( RAZ ) *S I N ( 6CLATT ) *C0S ( RL0NT ) -C0S ( RAZ ) *S I N ( RL0NT ) 

S I N ( RAZ ) *S I N ( GCLATT ) *S I N ( RL0NT ) +C0S ( RAZ ) *C0S ( RL0NT ) 
-SIN(RAZ)*C0S(GCLATT) 

-C0S( GCLATT )*C0S(RL0HT) 

-C0S^GCLATT)*SIN(RL0NT) 

-SIN(GCLAn) 



Figure A-2.- Analytic Drao Control Entry Guidance System block diagram. 





DRAG = XLFAC/^1.+XL0O2 


V<4572 


r>670*i.6( 



19126. 20- 1.625*V 



HS = 6248.4 


RGP(l) = RG(1)+0BAR 

RGP(2) = RG(2)-SIGN(RTURN,RG{2)) 

RGP(3) ° RG(3) 

DT * / AMAX 1 (0 . » RGP (1 ) Z+RGP ( 2 )2-RTURn2 ) 

DTH * TAN-1(RTURN/DT) 

DVHEAD = TArrl(-RGP(2)/(-RGP(l)-SIGN(DTH,RG{2))} 
U = I DVHEAD 1 


U » 2it-U 
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ARC = RTURN*'J 

RGP(l) = -DT*C0S(DVHEAD) 

RGP(2) = -DT*SIN(OVHEAD) 

RG{3) = RG(3)-RTE 

TARE(J) = RG(J)-RGP(J) ; J = 1,2,3 
UXYZE * [RG]/|RG! 

UTARE * ITAREI/ITARE! 

CTH = [UXYZEl* lUTAREl 
TRANGE = ATK*C0S-1(CTH)+ARC 
TA(J) = -RGP(J) ; J = 1,2,3 
[TA] = [RECll [TA] 


CIGARCl.l) = -SIN(GCLA1>C0S(RL0:O 
CI6AR(1,2) = -SIN{GCLAT>SIN(RL0N) 
CIGAR(1,3) = C0S(GCLAT) 

CIGAR(2,1) = -SIii(RL0N) 

CI6AR(2,2) = C0S(RL0:O 
CIGAR(2,3) = 0. 

CI6AR(3,1) = -C0S(GCLAT)*C0S{RL0N) 
CIGAR(3,2) = -Cg)S(GCLAT)*SIN(RL0N) 
CIGAR(3,3) = -SIN(GCLAT) 

[TAP] = [CIGAR] [TA] 


TAP(J)=SIGN{i:^°,TAP(J))j 

I 


PSIET * TArrl[TAP(2)/TAP(D] 
OELAZ * PSIE-PSIET 









Figure A-2. - Continued. 
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Figure A-2.- Continued. 















































D23 = -GSMVBB^/VSAT^ - 1. )/{1.4*C0S(BA)) 

AK = -3,5*D23/VBB 

C23 = -AK/3352.8 

C22 - 15849. 6 *C23 

C21 - 023 - C22*VBB - C23*VBb2 

nil = -.5'N((C21 + C22*VBB + C23*VBB^)/(C21 + C22*V + C23*V^))/C23 
SQ = 4.*C2o-C21 - C222 


SQQ = , 

H22 = C22*(TAN“^2.*C23*VBB/SQQ + C22/SQQ) 

-TAN“1(2.*C23*V/SQQ + C22/SQQ))/(C23*SQQ) 


SQQ = ASQ 

El = (2.*C23*VBB + C22 - SQQ)/(2.*C23*VBB + C22 + SQQ) 
E2 = (2.*C23*V + C22 - SQQ)/(2.*C23*V + C22 + SQQ) 
n22 = C22*LN(E1/E2)/(2.*C23*SQQ) 


RFF = nil + TT22 


/CG > V 


V < 6248.4 


RFF = 0. 


= .5*TEMP*LNJ;(VCg2 - VSAT^)/(VBB" 
= (VCG2 - VQ2)/(2.*ALFM) 


VSAT^)) 


S, T, 0; IMP 


Figure A-2.- Continued. 
















0 




DREFP = T2 

RDTREF = -2.*HS*DR£FP/V + C4 + C5*V 




Figure A-2.- Continued. 







Figure A-2.- Continued. 
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<^V < 2316.48 
N0V 


CAGI = l./(V^ + 2.*GS*HS) 

RDTREF = (-HS*(2.*DREFP*v - C11*V ) 


YES I C16 = 2.*C16 
n C17 = .9*C17 



C4*V^ + C5*V^)*CA6I 


I0F63 > 0 


V20LD * V0lO2 • 
RDC = (RDTREF - 
IDFG3 = 1 

152.4 

RDT0L2)/152.4 


U 


V > V20LO 


RDTREF » RDTREF - RDC*(V - .'20LD) 


Figure A-2.- Continued. 










Figure A-2.- Continued. 







Figure A-2.- Concluded. 
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XIII. APPENDIX B 
Digital Autopilot 

The Digital Autopilot (DAP) is designed to automatically fly the 
space shuttle orbiter from deorbit to the Icrminal Area Energy Manage- 
ment (TAEM) interface which occurs at an altitude of approximately 
21.3 km (70 000 ft) with a velocity of 457.2 m/sec (1500 tt/sec) . The 
DAP directs both the reaction control system (RCS) , and the aerodynamic 
control surfaces. 

The speed brake (6 sb) body flap (6gp) deflection schedules 
are shown in figure B-1, where figg is determined from a preset velocity 
schedule and 6gp is dependent on the center-of-gravity location. 

Figures B-2 through B-11 are block diagrams of the various elements of 
the DAP. Two types of signal limiting filters are used in this autopilot. 
The first type is illustrated below; 



This filter limits the value of the quantity (A-B) to be between x and 
y. The second type, called a "hysteresis filter," can appear in one of 
t io ways : 


b -a 
^ 



-b -a 


err 

m 

a b 

A » 

sa 

wm 



( 1 ) 


( 2 ) 







As A increases from zero, "err" remains zero until point "b" is 
reached. At this time, "err" becomes the value indicated (either a 
constant value if filter is type 1 or equal to A ii filter is t>'pe 2), 
and contin" > as A increases. As A starts to decrease, it remains 
the value indicated until point "a" is reached, where "err" becomes 
zero again. A similar situation would exist for an A decreasing 
from zero. 

TTie elevens are used for both elevator (dg) and aileron (6^) 
functions. The elevator command block diagram is shown in figure B-2. 

The aileron functions in one of two ways depending on the flight regime: 
for a < 18* and M < 5, the aileron is used for roll attitude (<}!) control 
figure B-3(a) ; when these conditions are not present, the ailerons are 
used for turn coordination, figure B-3(b). If the orbiter has a lateral 
center-of“gravity offset, the number of positive yaw (roll) thruster 
firings will not equal the number of negative yaw (roll) thruster firings 
due to the induced sideslip. By counting the number of positive and 
negative yaw and roll thruster firings, it is possible to establish the 
"steady state" aileron deflection required to offset this induced sideslip 
This is the role of the "up-down counter" shown in figure B-4. The number 
in parentheses in the block diagrams are the expressed values in English 
units. Figure B-S shows that the commanded left and right eleven 
deflections are functions of fig, and rudder (Sj.), 

figure B-6, is used for turn coordination when the aileron is used for 
roll control. If the ailerons are being used for turn coordination, the 
nidder is inoperative. 
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The pitch RCt)» figure 3-7* is operative for <1 less than 9?8 Pa 
(20 psf). In this regime it is used, along with the elevator, for 
longitudinal control. 

The roll RCS, figure B-8, is operative for q less than U79 Pa (10 
psf) and is used, along with the ailerons, for turn coordination. 

The yaw RCS, figure B-9, is oper«*ive throughout the entry until 
TAiM and serves one of two purposes impending on the flight conditions. 
If the ailerons are used for attitude control, the yaw RCS, 
figure B-9 (a), aids the rudder in maintaining turn coordination. 

If the I nditions are such that the ailerons are used for turn coor- 
dination, the yaw RCS, figure B-9(b), is used for roll attitude (^) 
control . 

lo integrate the linear first order differential equation in the 
control system, a convolution technique is used. his is a one-pass 
scheme that has demonstrated a high degree of accuracy in other real- 
time simulations, including piloted simulations. To illustrate, refer 
to figure B-10, which shows a typical first order system, x(t) W x(t) 


U(t) where U(t) is the forcing function. 


x(t) = e“W^ x(0) 




The solution is 
U(t) dx 


The convolution technique is a numerical method based on a Taylor 


series approximation (first two terms) of tie forcing function, C, am 
results in the following difference equation: 
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x(tj^+h) * P(h) x(tjj) + Q(h) ’-<(tj^) 
-Wh 

where P(h) * e 

«(h) • (h) , qj (»>1 


O(t^) = 


hr 


q^(h) » J dT » (1 - e“^)/W = (1 - P)/W 


qzCh) 


. Yte-"'" - ’> dT = - e-*^ . Whj/W^ • (h 


qi)/w 


0^ 


The control actuators, figure B— 11, are handled the sane way, except 
that provisions are nade for both position and rate limits. 

The RCS model uses the following equations to account for aerodynamic 
interference : 

Lrcs - hij HRJP - RJK) (YJP - YJri) 1 

«RCS ' t(PJP) - (PJN) + (YJP + YJN) 

+ (RJP + RJ!I) 

MpCS “ J(YJP - YJIJ) Kjj + (RJP - RJIl) Cj,J 

The values for the coefficients arvi shown in Table 3-1 . 
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SYMBOLS 


ir cine ter 

Uni t 

Oef i ni tion 

ay 

2 

m/sec 

side acceleration at center of gravity 

=e 

m 

elevon reference chord 


m 

rudder reference chord 

^he 

n,de 

elevon hinge moment coefficient 


n.d* 

^ (rudder hinge moment) / 3$ 

'“'r 

n.d. 

3 (rudder hinge moment) /36r 


n*d* 

rolling moment coefficient due to yaw RCS 

CmL 

n.d. 

pitching moment coefficient due to roll RCS 

Cmn 

n«d. 

pitching moment coefficient due to yaw RCS 

^NL 

n.de 

yawing moment coefficient due to roll RCS 

DEMX 

deg/sec 

maximum elevon rate 

DRMX 

deg/sec 

maximum rudder rate 

EP 

n.de 

pitch RCS error signal 

Er 

n*d. 

roll RCS error signal 


n«d« 

yaw RCS error signal 

f(6e) 

deg 

function of 6e used to limit 5_ _ 

a,c 

g 

m/sec^ 

acceleration. due to gravity 

h 

sec 

integration step size 

V 

N*m 

elevon hinge moment 

V 

N-m 

rudder hinge moment 


n.d. 

rolling moment RCS amplification factor 

%D 

n«d* 

pitching moment RCS amplification factor 
due to down firing Jets 

*Siu 

nrd. 

pitching moment RCS amplification factor 
due to up firing Jets 
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Kp 

n.d. 

di leron gain 

N* 

n*d« 

elevator gain 


n.d. 

rudder gain 

•■RCS 

N*m 

rolling moment due to RCS 

•“RJ 

N*m 

ideal rolling moment due to firing of 1 roll 
jet 

M 

n.d. 

Mach number 

Mpj 

N»m 

ideal pitching moment due to firing of 1 pitch 
jet 

MrCS 

M*m 

pitching moment due to RCS 

NrCS 

N«m 

yawing moment due to RCS 

Nyj 

N«m 

ideal yawing moment due to firing of 1 yaw 
jet 

p 

deg/sec 

roll rate 

p 

n.d. 

convolution ccef f i cient 

PJN 

n.d. 

number of negative pitch jets firing 

PJP 

n.d. 

number of positive pitch jets firing 

q 

deg/ sec 

pitch rate 

q 

Pa 

dynamic pressure 


sec 

convolution coefficient 

*»2 

sec^ 

convolution coefficient 

a 


vector of convolution coefficients 

r 

deg/sec 

yaw rate 

r' 

deg/sec 

r-(l80 g sinO cos 0 )/ttVj^ 

RJN 

n.d. 

number o. negative roll jets firing 

RJP 

n.d. 

number of positive roll jets firing 

8 


' Laplacian or<*rator 

Se 


eleven refer nee area 

Sr 


rudder reference area 
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t 

sec 

time 

tk 

sec 

time at the k sample 

U(t) 


convolution forcing function 

u 


vector of forcing function terms 

U(t) 


dU/dt 

V 

m/sec 

atmospheric relative velocity 

Vr 

m/sec 

earth relative velocity 

W 

sec“^ 

filter root 

x(t) 


convolution state vari»i ’e 

x(t) 


dy./dt 

y 

m 

lateral offset 

YJN 

n.d. 

numher of negative yaw Jets firing 

YJP 

n.d. 

number of positive yav Jets firing 

a 

deg 

angle of attack 

«c 

deg 

commanded angle of attack from guidance 
system 

6 

deg 

angle of sideslip 

e 

deg 

pitch angle 


deg 

roll angle 

♦c 

deg 

commanded roll angle to control system 


deg 

aileron deflection 

®a,c 


commanded aileron deflection 

6a,UD 

deg 

commanded aileron deflection from 
up-down counter 

6bF 

deg 

body flap deflection 

6e 

deg 

elevator deflection 

®e,c 

deg 

commanded elevator deflection 


deg 

left eleven panel deflection 

^ei,c 

deg 

command left elevon panel deflection 
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e,£m 

deg 

maximum change in elevon command 
allowed by rate limit 

'er 

deg 

right elovon panel deflection 

'<^r,c 

deg 

commanded right elevon panel 
deflection 

e,t 

deg 

Initial elevator setting 

'r 

deg 

rudder deflection 

r,c 

deg 

commanded rudder deflection 

r.llEi 

deg 

maximum change in rudder command 
allowed by rate limit 

3E 

deg 

speed brahe deflection 


sec 

variable of integration 
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q(Pa) 

0 

119.7 
239 .^ 

478.8 
718.2 
957.6 


TABLE B-I.- INTERFFRENCE RCS VALUES 


Jet Moment Value 






11185.5 





Mpj 



38325.6 





%J 



33878.8 




kl 

•^Mu 

kmd 

Kn 

cln 

Cmn 

CmL 

^NL 

746 

1.0 

.740 

1.02 

-.624 

0 

.130 

-.141 

688 

1.0 

.678 

1.02 

-.r3 

00 

0 

• 

.161 

-.115 

630 

1.0 

.616 

1.02 

-1.069 

.076 

.192 

-.111 

533 

1.0 

.541 

1.02 

-1.069 

.114 

.230 
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BODY FLAP, 6bf. SCHEDULE 



Figure B-1.- Booy flap and speed brake schedules. 





Figure Elevator conmand block diagram. 















Figure B-3.- Concluded. 
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Figure B>4.- Up-down counter block diagram. 













Figure B-6.- Rudder command block diagram. 
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Figure B-8.- Roll RCS error signal block diagram. 
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Figure B-9.- Yaw RCS error signal block diagram. 










Figure B-9. Concluded. 
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x(t|j + h) « P x(t|j) + Oi U(tk) + q2 U(t|^) 

x(t|t + h) * U(tk + h) - W x(t|^ + h) 

where 
P » e“^‘' 

q, = 0 - P)/W 
qg = (h - q^)/M 


Figure B-lO.- First order system. 















